Dynamic screening and plasmon spectrum in bilayer graphene 
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We have theoretically studied the collective response properties of the two-dimensional chiral elec- 
tron gas in bilayer graphene within the random phase approximation. The cooperation of external 
controlling factors like perpendicular electric bias, temperature, doping, and substrate background 
provides great freedom to manipulate the dynamic dielectric function and the low-energy plasmon 
dispersion of the system. Intriguing situations with potential application are systematically explored 
and discussed. Extra undamped plasmon modes might emerge under electric bias. They have almost 
zero group velocities and are easy to manipulate. 

PACS numbers: 71.10.-w,75.10.Lp,75.70.Ak,71.70.Gm 



INTRODUCTION 



Experimental breakthrough in isolation of high-quality 
few-layer graphene by exfoliation and epitaxial growth 
has led to intense experimental and theoretical interest 
in graphene materials Recently, much interest has 
been focused on the AB-Bernal stacked bilayer graphene 
(BLG) for fundamental physics and application poten- 
tial in nanotechnology [2h23|. With its own special na- 
ture, BLG inherits some characteristics of the monolayer 
graphene (MLG) carrying chiral Dirac fermions. Intrin- 
sic BLG is identified as a zero-gap semiconductor with 
quadratic band dispersion in low-energy regime instead 
of the linear band dispersion in MLG. In such a case, 
the density of state (DOS) in BLG at the two nonequiv- 
alent Dirac points, called the K and K 1 points, is a con- 
stant, in contract to vanishing DOS in MLG. Hence, BLG 
shares some similar features with two-dimensional elec- 
tron gas (2DEG). Specifically, an energy gap between the 
conduction and valence bands can be easily opened and 
tuned by introducing an electrostatic potential bias be- 
tween the two graphene layers The bias modifies 
the parabolic band structure, and increases the DOS at 
the top of the valance band and the bottom of conduction 
band. 

Previous studies have clearly evidenced that Coulomb 
interactions play a significant role in graphene [24], HH . 
Similarly, electron-electron interactions in BLG can also 
lead to exotic phenomena. The many-body effects are 
crucial to understanding the transport and optical prop- 
erties of the system. A particularly fruitful phenomenon 
is the dynamic screening. The frequency dependent 
screening determines the elementary quasiparticle spec- 
tra as well as the collective modes. At zero tempera- 
ture, the screening properties and the plasmon spectrum 
in BLG have been studied a naly tically in the two- and 
four-band approximations |7H10l|. For systems at finite 
temperature and under electric bias, numerical calcula- 
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tions have been employed to study the properties in the 
two-band approximation [ll|, dl ■ Because of the distin- 
guished energy band dispersions and chiralities, Coulomb 
screening properties and collective excitations in BLG 
exhibit significantly different behavior from the MLG 
and conventional 2DEGs. Experimentally the plasmon 
spectra and damping properties in graphene structures 
have been studied by methods such as electron-energy- 
loss spectroscopy (EELS) [H41i|. 

Surface plasmon modes have been used for the rapidly 
developing terahertz (THz) technology [13, HH and stim- 
ulated emission of plasmon in graphene has been pro- 
posed to be used as THz laser at room temperature 32j ■ 
One efficient way to increase the net plasmon gain is to 
decrease the group velocity of plasmon in the system. 
The plasmon dispersion in BLG is similar to that in 
MLG and might be easier to manipulate with the help 
of applied electric bias between the two layers. In this 
paper, based on Ref . [l2| . we numerically study the dy- 
namic screening properties of the Coulomb interaction 
in BLG systems within the random phase approxima- 
tion (RPA) and consider a fair general situation in which 
both the temperature and the bias voltage are finite. We 
calculate the dielectric function of BLG e(q, w) at arbi- 
trary wavevectors q and frequency lj. The zeros of the 
real part e r give the dispersion of the plasmon modes; 
the imaginary part indicates their damping properties 
to single particle excitations; the imaginary part of 1/e 
is related to their optical spectral weight. Beyond from 
Ref.[l2j], here we have explored systematically and in a 
high accuracy the screening properties and the plasmon 
spectra in a wide parameter space, and in some param- 
eter regimes have found two extra plasmon modes with 
almost zero group velocities. 



II. MODEL 

Since there are four inequivalent carbon atoms, the 
BLG system should be described by the four-band model 
which gives a hyperbolic dispersion [l8| . In the four-band 
model, a split-off band is located at A sa 0.4 eV above the 
lowest conduction band and a sombrero shape of band is 



2 



formed at the bottom of the conduction band in the ex- 
istence of electric bias. If the Fermi energy (ps 45.4 meV 
under U = 60 meV for n = 10 12 cm -2 ) plus the thermal 
energy (pa 26 meV at T = 300 K) is much lower than 
A, the effect from the split-off band should be negligible. 
For a bias potential U — 60 meV, the height of the som- 
brero is SU = 0.5U(1 - A/VU 2 + A 2 ) « 0.3 meV HO. 
In systems with carrier density n higher than the critical 
density n c = U 2 /{-Kh 2 v 2 F ) (0.27 x 10 12 cm" 2 for U = 60 
meV), the Fermi energy is above the sombrero and the 
topology of the Fermi surface is not affected by the som- 
brero. In this case, the effect of the sombrero appears as 
a DOS broadening (with a width of the height of the som- 
brero) at the bottom of the conduction band. Since one 
part of the thermal effect is also similar to a DOS broad- 
ening with a width of the thermal energy, the effect of 
the sombrero can also be neglected if the thermal energy 
is much higher than the sombrero height. Therefore, the 
low-energy properties of the system can be qualitatively 
well characterized by the two-band parabolic approxima- 
tion which is valid in the range of density n (0.5 to 2 
xlO 12 cm -2 ), temperature T (4.2 to 300 K), and bias 
potential U (up to 60 meV) most interested in this study. 
Within such a picture, a pair of chiral parabolic electron 
and hole bands touch each other at the Dirac (or the 
charge neutrality) point, and each band has a four-fold 
degeneracy arising from spin and valley degrees of free- 
dom. Note that the quadratic band dispersion in the two- 
band approximation deviates from the hyperbolic band 
dispersion predicted by the four-band model even in the 
intermediate energy regime. The corresponding results 
in some situations might differ quantitatively from those 
in real systems. For the sake of completeness, some zero- 
temperature results of systems with low carrier density, 
where the sombrero effect might not be negligible, are 
presented in the following. 

The low-energy effective Hamiltonian describing elec- 
trons of moderate energies in the K valley of biased BLG 
reads as Qjl 
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In the first term, the wavevector k = (k x ,k y ) 
k± = k x ±ik y is measured from the K point; the effective 
mass is m* = ^/(2v F ) ~ 0.035too with 7 the interlayer 
tunneling amplitude inherent in the BLG system, vf the 
graphene Fermi velocity, and too the free electron mass. 
The second term arises from the electrostatic potential 
bias U between the two graphene layers separated by a 
distance d = 3.35A. The eigenenergy and eigenwavefunc- 
tion of the above Hamiltonian read tl2l 



k, i.e., tan 9^ =k y /k x , k = |k|, and ak indicates the ratio 
of the kinetic energy to the potential bias with tan ak = 
H 2 k 2 /(m*U). For such an energy dispersion the DOS of 
the system is 
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where g u is a constant degeneracy factor. Here g v = 4 
comes from the degenerate two spins and two valleys at 
K and K'. Under finite bias U, the DOS of the BLG 
diverges on the edge of the energy gap E = \U/2\. It has 
been shown that the static dielectric constant at q = 
in intrinsic BLG is much larger than that in MLG due 
to the presence of a finite DOS at the K point [ll|, [T2I ]. 
The bias might further modify the system's properties in 
a large range of variety. 

The wavevector and frequency-dependent dielectric 
function e(q, oj) tells the response of the system to a weak 
external perturbation, and determines a variety of many- 
body related terms such as self-energy, carrier lifetime, 
and mobility, as well as characterization of other excita- 
tions. In the RPA, it is given as 



e(q,w) = 1 - v q n{q,u) 



(5) 



where v„ 



; 2 /(2e £b<7) (with the background dielectric 
constant e^) is the Fourier transformation of the bare 
Coulomb interaction and the electron-hole propagator is 
originated from the bare bubble diagram [121 ] 



n( q>w ) = 4 lffk A '(q) 



A,A',k 



(6) 



Here f(x) is the Fermi function and the vertex factor 
reads 

|^' A '(q)| 2 = |(k + q,AV^|k,A)| 2 

= - [1 + AA' COS Q!k COS Q!k+q 

+ AA' sin ak sin ak+ q cos(2#k — 2#k+ q )] (7) 

When q = and q = -2k, \g£ A '(q)| 2 = ±(1 + AA'). 
Similar to unbiased BLG, the interband vertical and 
back scatterings are both forbidden but the intraband 
back scattering is allowed in biased BLG. For intra- 
band scattering with k _L k+q, we have |<7j^ A (q)| 2 = 
rj[l + cos(ak + ak+q)], which becomes zero in unbiased 
BLG El. 



= e ik-r 



(2) 
(3) 



sin(^ + ±±Att) 

cos(^ + ±±kny 2e 

where A = ±1 denotes respectively the conduction and 
the valence band. Here #k is the azimuth of the vector 



III. RESULT AND DISCUSSIONS 

The electron-hole propagator is composed of intra- 
(A = A') and inter-band (A = —A') components. The 
intraband component is expected to be similar to that 
in conventional 2DEGs except the effect of chirality and 
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FIG. 1: The negative intraband electron-hole propagator in 
unit of No versus frequency uj in unbiased/biased BLG at 
zero/room temperature. The imaginary part for q = 0.06, 0.8, 
and 2.5/cf is plotted in panel (a), (b), and (c), respectively, 
and the corresponding real part in (d), (e), and (f). The 
electron density is n — 10 12 cm -2 with the Fermi wavevector 
fci? = 1.77 x I0 8 m -1 . The corresponding Fermi energy is 
E F = 34 meV at U = and E F = 45.4 meV at U = 60 meV. 
Solid curves are for T = 0, U = 0; dashed T = 0, U = 60 
meV; dotted T = 300 K, U = 0; dash-dotted T = 300 K, 
U = 60 meV. 

deformation of energy band. The interband one, which 
can be manipulated by the bias voltage, modifies qual- 
itatively the screening properties of the system. Since 
the approximate energy spectrum we use in this study is 
isotropic, the obtained properties are also isotropic and 
we use q = |q| and k = |k| in the following discussion. 

In Fig. [TJ we plot the negative intraband propaga- 
tor in unit of iV = 2m* /tt (the DOS of intrinsic BLG) 
versus frequency uj for three typical wavevector q val- 
ues. Similar to the 2DEG result at zero temperature, 
the imaginary part is nonzero in the single particle con- 
tinuum < uj < uj u for < q < 2kp and oji < uj < lu u for 
q > 2k F , wither = \El p _ q -EU andw u = E^-E 1 ^. 
The derivative of the imaginary part is not continuous at 
the continuum edges and for < q < 2kp is also not con- 
tinuous at uj = uji besides uj = and uj u . When the bias 
U increases, the structures of the curves shift to lower 
energy because the energy band is narrowed and the 
group velocity of electrons at the bottom of the conduc- 
tion band decreases. When the temperature increases, 
the sharp edges become smoother and the nonzero range 
gets wider as expected. 

Different from the 2DEG result, as illustrated in 
FigHJc) for q > kp at U = and T = 0, the imagi- 
nary part has a sharp dip with a derivative discontinuity 
at uj m — q 2 /2m* between the edges due to the chiral 
nature of the wave functions @. The dip and disconti- 
nuity persist at finite temperature but become softened 
and disappear under finite bias voltage. 

The real part, presented in the lower panels of FigJlJ 
shows a sharp peak at uji and a sharp dip at uj u with an 



FIG. 2: The negative interband electron-hole propagator in 
unit of iVo versus frequency uj is plotted for the same param- 
eters set and arrangement as in Fig[T] 



extra peak near uj m for q > kp and moves in a similar 
way as the imaginary part with bias and temperature. 
However, the structure near uj m for both the imaginary 
and real part develops in a different way from those at 
oji and uj u when U and T increase. It remains at high 
temperature and removes off with the bias voltage while 
those at uji and uj u decay with temperature. On average, 
the variation range of the propagator decreases with q. 

In Figl^l we present the negative interband propaga- 
tor in unit of Nq versus uj for the same three q values as 
in Fig[TJ Its imaginary part and the single particle con- 
tinuum have a minimal energy limit uj\ = E\ p — E^_ 
at zero temperature or uji — 2E^ 2 at finite tempera- 
ture. For small q [Figj^a)] the imaginary part increases 
sharply and reaches a peak before decreases in a way 
~ 1/uj. For mediate q = 0.8/cf [FigHJb)] the main 
peak becomes round and smooth while a sharp peak near 
ujs = E\ p — E^ p+q appears under bias as shown in the 
inset of FigJSJb). This latter peak grows and becomes 
more visible as U increases. For q > kp a peak with 
discontinuity appears near uu m — q 2 /2m* at zero tem- 
perature in both biased and unbiased BLG Q as shown 
in FigHc). 

Corresponding to the continuum edges of the imagi- 
nary part, on the curve of the real part as shown in the 
lower panels of Fig[21 there is a peak near wi (uJ2) at low 
(high) temperature when the electron system is degen- 
erate (nondegenerate) . The peak near uji may become 
very sharp for small q in biased BLG because the bot- 
tom (top) of the conduction (valence) band becomes flat 
and the DOS diverges on the edge of energy gap. The 
overall contribution of interband excitation to the prop- 
agator increases with q and in a way ~ q 2 at small q. 

As we know, the electron-hole propagator reflects the 
electric polarizability of a many-body system screening 
a Coulomb potential. After being reduced by the back- 
ground dielectric constant £(,, it determines the dielectric 
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FIG. 4: Plasmon spectrum of intrinsic BLG (a) in vacuum 
with Eb = 1 at various temperature and (b) with various 
dielectric constant at room temperature T = 300 K. The 
shadow shows the electron-hole single-particle continuum at 
zero temperature. 



FIG. 3: The real part of the dielectric function e r vs fre- 
quency uj when et — 1 in variety of bias U, temperature T, 
electron densityn (xlO 12 m -2 ), and wavevector q (xlO 8 m _1 ) 
is plotted in (a)-(d), respectively. Parameters U = 30 meV, 
n = 10 12 m" 2 , q = 0.1 X 10 s m" 1 , T = 300 K are used if not 
specified in the panels. 



function e(q,uj) = e r + ia as indicated in Eq.Q. The 
zero of the real part e r gives the collective excitation of 
the system in the absence of external electromagnetic 
field. The imaginary part Ci gives the spectrum of the 
single particle excitation. In the presence of only intra- 
band single particle excitation as in conventional 2DEG, 
there exist maximally two plasmon modes, one acoustic 
mode of frequency uja within the single particle excita- 
tion and another optical mode u>o, because — Re[n(g,u;)] 
has only one dip below zero as shown in Fig[TJ The acous- 
tic mode is thus always overdamped with little spectral 
weight and not experimentally relevant. The contribu- 
tion from interband excitation introduces fine structures 
to e r near zero and at least two extra modes, Wg and io\ 
may emerge. 

In Fig(3J we display e r versus uj for various temper- 
ature T (a), bias voltage U (b), electron density n (c), 
and wavevector q (d). In the high frequency limit, the 
effect of polarization vanishes and £oo = 1. For intrinsic 
BLG where n = 0, the intraband polarization is negligi- 
ble and e r > at T — 0. There is no collective mode. In 
other cases, the e r versus w curve has a deep dip at uj u 
and a peak at u}\ or/and uj2- The competition in Eq.([5]) 
among value one, the intra-, and inter-band contributions 
to the polarizability determine its fine features. For typi- 
cal parameters, as illustrated in Fig[T]and[21 the effect of 
the intra- (inter-) band polarization decreases (increases) 
with q so we expect that e r mainly shows intra- (inter-) 
band characteristics for long (short) wavelength. In this 
study we are mostly interested in the screening and col- 
lective excitation properties of long wavelength in the sys- 
tem, and in FigJ3]we present the details only for q -C kp 
near e r = 0. 

The result for a system of U — 30 meV, n = 10 12 



cm~ 2 , and q = 0.1 X 10 8 m _1 is exhibited in FigJ^a). 
At low temperature T = 4 K (solid) when the system 
is degenerate, the intraband contribution dominates and 
a peak appear near uj\ — E\ f — E^_ q ~ 1E\ p = 74.4 
meV which is out of the panel. When the temperature 
increases to 200 K (dashed) and the system becomes non- 
degenerate, a sharp peak emerge near UJ2 = 2£ , ^ 2 — 30 
meV. Since this interband peak is located just below the 
optical plasmon energy ujo, it can introduce two extra 
roots and ix>\ to the equation e r = 0, i.e., two extra 
plasmon modes in the system. The extra plasmon mode 
of lower energy is out of the interband single particle 
continuum and can be only weakly damped. When the 
temperature increases further (dotted and dash-dotted) 
the enhanced intraband contribution shifts the peak be- 
low zero and the extra plasmon modes disappear. In the 
same time, loo increases with the temperature as more 
electrons (holes) exist in the conduction (valence) band. 

As shown in Fig(3Jb), the bias voltage can sensitively 
shift the interband peak and control the emergence of 
the extra plasmon modes. These modes have energies 
(frequencies) proportional to the bias voltage and group 
velocities close to zero. The increase of the electron den- 
sity enhances the degeneracy of the system and reduces 
the amplitude of the interband peak at 0J2 — 2E 1 ,„ as 

qj z 

plotted in FigJ3{c). The enhanced intraband contribu- 
tion at higher density also increases the frequency of the 
optical plasmon mode ujo- 

In Fig[3{d) we illustrate how the e r -u> curve develops 
with q at room temperature T = 300 K. Overall the 
intra- (inter-) band contribution to e simply decreases 
(increases) with q as shown previously in Figs[T] and [21 
but their effect on the plasmon spectrum is more compli- 
cated due to their competition with each other. When q 
increases, the intraband introduced e r dip becomes wide 
which usually results in the increase of the plasmon fre- 
quency. The interband peak of e r is located at the fixed 
energy ui2 — 30 meV and its amplitude increases with 
q. As a result, only when 0J2 is close to wo, the inter- 
band contribution can affect significantly the plasmon 
spectrum. 
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FIG. 5: Plasmon spectrum at zero temperature of (a) unbi- 
ased (b) U=60meV biased BLG in vacuum for density n — 0.1 
(solid), 0.5 (dashed), 1.0 (dotted), and 2.0 (dash-dotted) 
xlO 12 cm -2 . The single-particle continuum at n = 10 12 cm~ 2 
is also present (light/dark shadow for the inter-/intra-band 
part). 



The plasmon spectrum in intrinsic BLG (unbiased and 
undopcd) depicted in FigJH at various temperatures T 
(a) and for various background dielectric constant Et (b). 
At zero temperature, there is no carrier and no plas- 
mon mode in the system. At finite temperature, elec- 
trons (holes) are excited in the conduction (valence) band 
and two plasmon modes emerge. In the long-wavelength 
limit, their frequencies are proportional to y/T at high 
temperature. Similar to the 2DEG result, we also ob- 
serve a dispersion uo oc ^fq and oja oc q. The background 
dielectric constant can also be employed to modify the 
plasmon frequency as shown in Fig|4]Jb). The acoustic 
mode is not sensitive with £h but the frequency of optical 
mode decreases quickly with e^. 

In Fig. [SI the plasmon spectrum for various electron 
densities at zero temperature is illustrated in unbiased 
BLG (a) and in biased one with U = 60 meV (b). In the 
long- wave limit, q ~ 0, the dielectric function is dom- 
inated by the intraband contribution. The properties 
of the system at zero temperature is mainly determined 
by the group velocity of electrons at the Fermi energy. 
The plasmon spectrum of unbiased system is similar to 
that of 2DEG. If the electron density is not high, e.g. 
n < 10 12 cm~ 2 for U = 60 meV as shown in Fig(5jb), 
the Fermi group velocity of electrons decreases when the 
external electric field is turned on, due to the band defor- 
mation, and the plasmon mode becomes softened. With 
increasing q the interband contribution becomes more im- 
portant, which reduces the group velocity of the optical 
plasmon mode. In some cases, the group velocity of plas- 
mon can be close to zero, a favorite situation for the 
stimulated plasmon emission [32| . For large q, when the 
optical plasmon branch enters the interband single par- 
ticle continuum, i.e. uoo > wi, the effect of interband 
contribution decreases and the plasmon spectrum has a 
long tail in unbiased system or when the carrier density 
is high. However, in biased systems with low carrier den- 
sity, the effect of the interband contribution can be sig- 



nificant due to the flat band and the plasmon spectrum 
ends near where the intra- and inter-band single particle 
continua meet at uj u = uj\. 

The competition between the intra- and inter-band 
contributions may result in two extra plasmon modes for 
proper U at finite temperature when thermal excitation 
becomes important as previously discussed in Figj3] This 
is an interesting phenomenon because those modes have 
unique properties and might be used in nanotechnology. 
In FigEl we plot the plasmon spectra of a biased system 
with U = 30 meV at T = 300 K and n = 10 12 cm~ 2 for 
background dielectric constant = 1 (solid), 5 (dash- 
dotted), and 15 (dotted). At room temperature T = 300 
K, the electronic system is not degenerate and the plas- 
mon spectrum is in general similar to the one of intrinsic 
BLG as shown in Figrjjb). However, as shown in FigJ^] 
and [3l the interband contribution adds a sharp peak at 
u>2 which is around U for small q. As a result, the plas- 
mon spectra are deformed at ui2 and bifurcate in some 
cases where two extra plasmon modes Wg and uj\ emerge 
as shown in the inset of FigJBl 




0.0 0.4 0.8 1.2 1.6 

q (10 8 m 1 ) 

FIG. 6: Plasmon spectrum of biased and doped BLG with 
U = 30 meV and n = 10 12 cm -2 at room temperature 
T = 300 K in environment of different background dielec- 
tric constants £b- The inset shows the zoomed spectrum of 
the plasmon branches uio,^z, uj\, and uja for Sb = 1 at small 
q. The light (dark) shadow shows the interband (intraband) 
single-particle continuum at zero temperature. 

The emerged two plasmon modes have almost zero 
group velocities and their frequencies are proportional to 
the bias voltage. In addition, their energies are in the gap 
of the zero-temperature single particle continuum and the 
lower one is below the lower limit of single particle exci- 
tation, u>2, at high temperature. This suggests that the 
modes are undamped or weakly damped and have long 
lifetime. In FigJT] we plot the negative imaginary part 
of the dielectric function, which indicates the spectral 
weight of the collective modes, for different q. We see 
that the spectral weigh shows a negligible value (wide 
peak) near the energy of the acoustic (optical) mode uja 
(wq), indicating that the uja mode is damped and the ujq 
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FIG. 7: The spectral weight function versus frequency at var- 



ious wavevectors q (xlO 8 
FigH 



m ) for the system with ej = 1 in 



is weakly damped. In the same time, there is a sharp and 
high peak at uj i2 U for small q indicating that the W3 
and W4 modes are almost undamped. These undamped 
plasmon modes are similar to those in MLG but have 
lower group velocity. Their energies can be easily manip- 



ulated by the bias voltage. These undamped modes with 
almost zero group velocities then might be used in THz 
technology [32]. 

In conclusion, we have studied systematically the 
many-body response of electrons to external Coulomb 
perturbation and the plasmon spectra in a biased BLG. 
The vertical voltage bias opens a gap between the con- 
duction and valence bands and increases the DOSs at the 
band edges. As a result, the bias modifies the dielectric 
function greatly. In the long-wave limit, a sharp and con- 
trollable dielectric peak might appear at the energy equal 
to the band gap at high temperature when the system is 
nondegenerate or at the energy of double the Fermi en- 
ergy at low temperature when the system is degenerate. 
In some cases, two extra undamped plasmon modes ap- 
pear at energies close to the band gap energy and have 
almost zero group velocities. 
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